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Abstract 

We give an elementary proof of Sarason's solvability criterion for the 
Nevanlinna-Pick problem with boundary interpolation nodes and boundary 
target values. We also give a concrete parametrization of all solutions of such 
a problem. The proofs are based on a reduction method due to Julia and 
Nevanlinna. Reduction of functions corresponds to Schur complementation 
of the corresponding Pick matrices. 

1 Introduction 

In this paper we present an elementary solution of a version of the Nevanlinna- 
Pick interpolation problem with boundary interpolation nodes and boundary tar- 
get values. Such problems have been studied by many authors, beginning with R. 
Nevanlinna [H] in 1922, and they continue to attract interest, not only as natural 
questions in function theory, but also because they have applications to engineer- 
ing, particularly electrical networks and control theory (see [TJ page 812] for some 
references). The problem is to construct functions in either the Schur class S or the 
Pick class V (denned in Section [2]) subject to finitely many interpolation conditions 
on the values of the function and its first derivative at points on the boundary of 
its domain of analyticity. 

Nevanlinna, in his original paper [11] on the problem, gave a simple recursive 
technique based on a theorem of G. Julia to determine whether a boundary inter- 
polation problem is solvable, and if so, to describe all solutions. The method is 
analogous to "Schur reduction" for the standard Nevanlinna-Pick problem. Nevan- 
linna did not, however, give a criterion for the existence of a solution. In 1998 D. 
Sarason, after referring to subsequent papers by M. G. Krein, J. A. Ball, J. W. 
Helton, A. Kheifets and several others, wrote [12] "It is unclear, at least to this 
author, to what extent the preceding treatments can be extended so as to clarify 
Problems dNPS and dNPC, for example, and to find a criterion for the solution 
set to be non-empty." He went on to give a succinct criterion, after converting the 
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interpolation problems in question to moment problems on the circle. He also gave 
a description of the solution set in terms of measures on the unit circle that solve 
a form of truncated trigonometric moment problem. 

The elementary reduction technique of Julia and Nevanlinna has the virtue that 
it assumes only a good first course in complex analysis. It does not require any 
operator theory or even the notion of Hilbert space, and so may be particularly 
suitable for working engineers. In this paper we show that the technique can indeed 
be used to prove Sarason's solvability criterion, and furthermore to parametrize all 
solutions of Problems dNPV and dNPS (Sarason's terminology, explained below) 
in a straightforward and concrete way. Nevanlinna's own parametrization (TTJ Satz 
I] is marred by an oversight, which we correct in Theorem 12.21 We also remove 
an unnecessary restriction (distinctness of the interpolation values) and we incor- 
porate the condition for indeterminateness. It transpires that the parametrization 
problem is more delicate and even more interesting than Nevanlinna realized. 

There is a substantial literature on boundary interpolation. Besides the papers 
discussed above we mention particularly papers of J. A. Ball and J. W. Helton [3], 
D. R. Georgijevic [7J and V. Bolotnikov and A. Kheifets [5j, and the books of J. 
A. Ball, I. C. Gohberg and L. Rodman [2], and of V. Bolotnikov and H. Dym [I]. 
These authors variously make use of Krein spaces, moment theory, measure theory, 
reproducing kernel theory, realization theory and de Branges space theory. They 
obtain far-reaching results, including generalizations to matrix-valued functions 
and to functions allowed to have a limited number of poles in a disc or half plane; 
see [2j [5] both for results and for many references to the literature. We shall 
comment briefly throughout on the points of contact of their work with this paper. 
An elementary approach also has its virtues, and provides different insights into 
the subtleties of boundary interpolation. 

The paper is organised as follows. In Section [2] we state the problem and two 
theorems: Sarason's solvability criterion and the corrected version of Nevanlinna's 
parametrization of all solutions. We also briefly describe some other authors' 
parametrizations. In Section [3] we present an extension of the Julia-Nevanlinna 
reduction method and the inverse process of augmentation. We prove the key fact 
that Nevanlinna reduction and augmentation preserve the Pick class. In Section 
H] we give a simple example which both illustrates Nevanlinna's procedure and 
provides a counterexample to [ITJ Satz I]. In Section [5] we prove that Nevanlinna 
reduction of functions corresponds to Schur complementation of Pick matrices 
(Theorem 15. 3p . In Section Owe give the promised elementary proof of Sarason's 
result, Theorem 12. 11 In Section [7] we prove the parametrization theorem, Theorem 
I2.2[ and show that it can be interpreted as a linear fractional parametrization or 
as a continued fraction expansion (Corollaries 17.21 and 17. 3j) . In Section [8] we sketch 
the corresponding results for boundary interpolation in the Schur class. 

We are grateful to Joseph Ball and Vladimir Bolotnikov for drawing our atten- 
tion to several relevant papers. 
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2 The interpolation problem and two theorems 

We denote by II the upper half plane {z e C : Im z > 0}, and by D the open 
unit disc {z : \z\ < 1}. The Pick class V is the set of functions analytic and 
with non- negative imaginary part in II. The Schur class S is the set of functions 
analytic and bounded by 1 in D. We study interpolation problems for functions in 
V and S. Let us state the version for V . 

Problem dNPV: Given distinct points X\, . . . , x n £ R and numbers wi, . . . , w n £ 
R and v±, . . . , v n > 0, determine whether there is a function f £ V satisfying 

f(xj) = Wj, f'(xj) = vj for j = 1, . . . , n, (2.1) 

and if there is, describe the class of all such functions f . 

The Pick matrix for this problem is the matrix M = [?riy]^- =1 where 




if i = 3 
if i ± j- 



For simplicity we shall understand the conditions (12. ip . as Nevanlinna did, in the 
straightforward sense that / is analytic at each interpolation node, though Sarason 
interpreted them in a weaker sense. The simplification is justified by the fact that 
Problem dNPV has a solution in the weak sense (values in conditions (12. ip are 
non-tangential limits) if and only if it has a solution in the strongest possible sense 
(the function is rational and analytic at the interpolation nodes) - see Note 16.61 
and Proposition 16.71 

Sarason [12] studied the equivalent problems for S and the class C of analytic 
functions with non-negative real part in D. We state the problem dNPS in Section 
[8] below. 

A first guess, in analogy with Pick's theorem for interpolation at interior points 
of © or II, might be that there exists a solution of Problem dNPV if and only 
if the Pick matrix M for the data is positive. However, there is a simple reason 
why this cannot be so - the fact that a non-constant function in V has a positive 
derivative at any point in R at which it takes a real value (Proposition 13.11 (1)). 
Consider the following problem dNPV: 

Find f £ V such that f(xi) = f(x 2 ) = and f'(xi) = 0, /'(x 2 ) = 1> where 
Xi ^ x 2 E R. 

The constraint f'{x\) = forces / to be constant, and so there is no such /, 
but the Pick matrix for these data is 



which is positive. One might argue that the condition f'(xi) = is unnatural, and 
hope that positivity of M would suffice in the case of strictly positive Vj. However, 
it is not so: the above example can be modified to produce a counter-example with 
three interpolation nodes and all Vj > 0. 
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It is an intriguing feature of boundary interpolation problems that positivity 
of the Pick matrix is necessary but not sufficient for solvability. As far as we know 
this was first pointed out by B. Abrahamse [TJ page 812]. 

Here is Sarason's solvability criterion (but for the replacement of C by V) . We 
say that / is a solution of Problem dNPV if / G V, f is analytic at x%, . . . , x n and 
/ satisfies the conditions (12.11) . We say that Problem dNPV is solvable if it has at 
least one solution, is determinate if it has exactly one solution, and is indeterminate 
if it has at least two solutions. The problem has a convex solution set, and so has 
zero, one or infinitely many solutions. A positive semidefinite matrix is said to be 
minimally positive if it majorises no non-zero positive semidefinite diagonal matrix. 

Theorem 2.1 A problem dNPV is solvable if and only if the corresponding Pick 
matrix is positive definite or minimally positive. 

The solution set of the problem, if non-empty, contains a real rational function 
of degree no greater than the rank of the Pick matrix. 

The problem is determinate if and only if the Pick matrix is minimally positive. 

Our proof of Sarason's Theorem will be by means of an induction, which also 
leads to the following concrete parametrization of solutions. 

Theorem 2.2 Suppose that Problem dNPV has a positive definite Pick matrix. 
There exist 

Si, . . . ,s n G RU {oo}, t 1 ,...,t n >0, Ui,...,y n GKU{oo} 

such that the general solution of Problem dNPV is f = f\ where the functions 
f n +i, ■ ■ ■ , fx are given recursively by: 

(1) fn+i is any function in V that is meromorphic at Xi,...,x n and satisfies 
fn+i(xj) ^Vj, 3 = l,-..,n; 

(2) fk is the augmentation of fk+i at x^ by s/., tj, for k = n, . . . , 1. 
The quantities Sj,tj, yj for j = 1, . . . , n are given explicitly by 

■v o //, //.;' • 1 (2.2) 

where Wj, Vj for 1 < k < j < n and for 1 < j < k < n + l are given by equations 
flUD to flUD below. 

The notion of the augmentation of a function at a point of the real axis will 
be defined in the next section, and the equations (17.11) to (17.41) that define the 
quantities Sj, tj, yj are simple piecewise rational expressions. In this way we obtain 
a parametrization of all solutions of Problem dNPV as a linear fractional transform 
of a function f n+ i in the Pick class, which however is not quite free: it has to be 
meromorphic at the interpolation nodes and it must satisfy a single inequation at 
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each node. The linear fractional form of the parametrization is stated in Corollary 
17.21 It can also be stated in continued fraction form: see Corollary 17.31 

Theorem 12.21 is a corrected version of [HI Satz 1.2], where the inequations 
f n+ i(xj) ^ jjj were omitted on account of an oversight discussed further in Note 

El 

A more general (but accordingly less explicit) result is proved in [21 Theorems 
21.1.2 and Corollary 21.4.2] with the aid of realization theory. In this version one 
seeks a rational matrix-valued function / with non-negative real part and with at 
most k poles in II satisfying matricial versions of the interpolation conditions ( 12. ip . 
On specialization to scalar functions we obtain the following. If the Pick matrix 
of Problem dNPV is invertible and has at most k negative eigenvalues then the 
solution set consists of all functions of the form 

}{z) = (Qn(z)g(z) + Q 12 (z))(Q 21 (z)g(z) + Q 22 {z)Y 1 

where the 0jj are given by an explicit formula in terms of the data, and g ranges 
over the scalar rational functions that are bounded by 1 on II and are such that 
®2i<7 + @22 has simple poles at the interpolation points x\, . . . ,x n . The italicised 
condition is equivalent to a set of inequations g(xj) ^ i/j, as in Theorem 12.21 

Another parametrization is given in [51 Theorem 4.3] within the framework of 
the Ukrainian school's Abstract Interpolation Problem. The authors allow finitely 
many derivatives to be prescribed at each interpolation point. The proofs depend 
on operator theory in de Branges-Rovnyak spaces. 

J. A. Ball and J. W. Helton have developed a far-reaching theory of interpola- 
tion which they call the "Grassmanian" approach; it makes use of the geometry of 
Krein spaces. It gives a unified treatment of many classical interpolation problems, 
including the ones studied in this paper - see [31 Theorem 6.3]. Inevitably, one 
pays for the generality with a less concrete parametrization. 

Some terminology: for a matrix M, "positive" means the same as "positive 
semi-definite", and is written M > 0, whereas M > means that M is positive 
definite. In a block matrix 

" A B 
C D 



M 



where A is non-singular, the Schur complement of A is defined to be D 
By virtue of the identity 



CA~ X B. 



M 



' A 


B ' 




1 


" 




C 


D 




CA~ l 


1 





A 

D - CA- X B 



1 A- l B 
1 



(2.3) 



it is clear that rank {D — CA l B) = rank M — rank A. Furthermore, if A > 
then M > if and only if C = B* and D - CA~ l B > 0. 



3 Reduction and augmentation in the Pick class 

We shall need the following basic properties of the Pick class V. 
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Proposition 3.1 Let f be a non- constant function in the Pick class and letx E R. 



(1) If f is analytic and real-valued at x then f'(x) > 0. 

(2) If f is meromorphic and has a pole at x then f has a simple pole at x, with 
a negative residue. 

Proof. (1) Suppose /(£ + i^) = u + iv with £,rj,u,v real: since v > on II 
and v(x) =0 we have v v {x) > and hence, by the Cauchy-Riemann equations, 
u i{x) > 0. Furthermore the restriction of v to a neighbourhood of x in R attains 
its minimum at x, and so v% — at x. Hence f'(x) = (u^ + iv^)(x) > 0. 

Suppose f\x) = 0: we can write f(z) = f(x) + (z — x) u g(z) in a neighbourhood 
of x, with g analytic, g{x) ^ and v > 2. Pick a neighbourhood of x on which g 
is bounded away from and arg g lies in a small interval J about arg g(x). For 
small enough r > and < # < 7r, if^ = x + re 10 then 

arg (/(z) — f(x)) = u6 + arg g(z) E vQ + I. 

As 6* goes from to 7r the interval u6 + I rotates at least once around the circle, 
and so Im f(z) < at some point in II, contrary to the fact that f E V. Hence 
f'{x) > 0. 

(2) Observe that if / is a non-constant function in V then —1/fEV. If / has a 
pole of order v at x and leading term R/z u then —1/f has a zero of order v at x 
and leading term —z u /R. By (1) i/ = 1 and — 1/R > 0. □ 
Nevanlinna used a recursive technique for eliminating interpolation conditions 
at points on the real line. The reduction procedure (in the case of a function 
analytic at an interpolation point) is due to G. Julia [S|. 

Definition 3.2 (1) For any non-constant function f E V and any x E R such 
that / is analytic at x and f(x) E R we define the reduction of f at x to be the 
function g on n given by the equation 

1 1 

g[Z) ~~f{z)-f{x) + f'{x){z-x) (3 ' 1} 

For f E V having a pole at x E R with residue R we define the reduction of / at 
x to be the function g on n given by 

g(z) = f(z)-R/{z-x). (3.2) 

(2) For any g E V, any x E R such that g is meromorphic at x and any a E 
RU {oo}, ax > 0, we define the augmentation of g at x by ao, a\ to be the function 
/ on n given by 

1 1 

g(z) if a E R 



f(z)-a ax(z-x) ^ 

f{z) = g(z) r if a = oo. 

aAz — x) 
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The former of the equations in (13.31) can of course be written 

f(z) =a + 1 1 (3.4) 

Note 3.3 If / is a non-constant [real] rational function in V of degree N and / is 
analytic at x G R then its reduction at x is a [real] rational function in V of degree 
at most N — 1 . 

If g G is [real] rational of degree N then the augmentation of g at any point 
of analyticity in R by ao G R, «i > is [real] rational of degree N + 1, and the 
augmentation of g at any pole in R is [real] rational of degree N. 

The important property of the operations of reduction and augmentation is 
that they preserve the Pick class, as we now show. 

Theorem 3.4 Let x G R. 

(1) If a non-constant function f G V either has a pole at x or is analytic and 
real-valued at x then the reduction of f at x also belongs to V and is analytic 
at x. 

(2) If g G V is meromorphic at x and a G R, a,\ > then the augmentation f of 
g at x by ao, ai belongs to V, is analytic at x and satisfies f(x) = ao, f'{x) < 
a\. In fact if g is analytic at x then f'(x) = a\, while if g has a pole at x of 
residue R then f'(x) = ai/(l — a±R) and hence f'{x) < a\. 

(3) If g G V is meromorphic at x and a\ > then the augmentation f of g at x 
by oo, ai belongs to V and has a pole at x with residue — 1/aj. 

Julia [8] stated and proved necessity in (1) in the case that / is analytic at x. 
Nevanlinna [11] gave a different proof. E. Wigner [14], over 30 years later, stated 
the analogous result for the smaller class of U R functions" , that is, functions in the 
Pick class that are meromorphic on the entire plane and real on the real axis. He 
gave a more elementary proof, which he attributed to Schiffer and Bargmann. Some 
subsequent authors have followed Wigner in calling (1) the "Schiffer-Bargmann 
lemma" . 

The use of augmentation at a pole will be important in our treatment. Nevan- 
linna seems not to have considered it; we discuss the consequence in Note 16.41 
below and in Section [71 The fact that, in boundary Nevanlinna-Pick interpolation, 
positivity of the Pick matrix is insufficient for solvability is well explained in terms 
of the strict inequality u f'(x) < ai" for augmentation at a pole. 
Proof. (1) Suppose (Case (i)) that / has a pole at x with residue R. By 
Proposition 13.11 R < 0. The reduction g of /, given by equation (13.21) . is analytic 
at x; hence there is a neighbourhood U of x such that Im g is harmonic and 
continuous on II U U. 
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Let e > 0. On U n R we have Im g = Im / > 0, and so by continuity there 
is a neighbourhood V(x) of x on which Im g > —e. On the other hand, for real t 
distinct from x (including t = oo) there is a neighbourhood V(t) of t in C U {oo} 
such that for z G V{t) 

R 



Im I | > -e 

z — x 



and hence, for z G V(£) R II, 



Im (7(2:) = Im /(z) — Im > — e. 

Julia gives a geometric argument to show that g £ V, whereas Nevanlinna con- 
cludes: it follows by a well-known application of the Poisson integral that Im g > 
on n. 

Here is an alternative, more elementary, ending to the argument based on the 
idea of Schiffer and Bargmann p3l II. 8a]. 

Suppose that g V, so that there exists z' G II such that 

p = f -Im g(z') > 0. 

Choose e = \p. By compactness K U {00} is covered by finitely many of the 
V(t), and hence there exists 5 > such that Im g(z) > —\p for all z such that 
< Im z < 5. Define a contour C = C\ U C 2 in II as follows. Let r > 2\R\/p. We 
take C\ to be the portion of the circle about x of radius r that lies in the half-plane 
{Im z > and C2 to be the segment of the line {Im z = |5} lying between the 
two intersections of the line with C\. We can assume r so large and 5 so small 
that z' lies inside C . Then we have — Im g < |p on C, but — Im g(V) = p. This 
contradicts the maximum principle for harmonic functions, and so g G V . 

Now consider Case (ii): / is analytic and real- valued at x. Here f'{x) > and 
the reduction g of / at x is given by equation (13. ip . The function — l/(f — f(x)) is 
in the Pick class and has a simple pole with residue —l/f'(x) at x. Case (i) shows 
that g G V . 

(2) When / is the augmentation of g at x by a G R, ai > we have 



1 



0V* 



f(z)-a ai(z-x)' 

The right hand side is in the Pick class and so in turn are the left hand side, / — 
and /. Furthermore, from equation ([31 



/(*) = «o + n 7 - w : (3-5) 

Since g is meromorphic at x and is in V, g is either analytic at x or has a simple 
pole at x with negative residue, by Proposition 13.11 In either case —{z — x)g(z) 
is analytic and non-negative-valued at x, and hence the denominator of the right 
hand side of equation (13.51) is analytic and (since af 1 > 0) strictly positive at x. 
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Hence f(x) = do- If is analytic at x then f'(x) = a±, while if g has a pole of 
residue R at x then the denominator in equation (13. 5p takes the value a^ 1 — R at 
x, and so f'(x) = a\j (1 — d\R). Since a\ > and i? < we have 1 — a\R > 1, and 
so ai/(l — aii?) < ai. 

(3) The statement is immediate from the second equation in (13.31) and the fact 
that the function —l/(z — x) is in the Pick class. □ 

Note 3.5 A function g obtained by reduction is not a general element ofV. 

For if g is the reduction of a Pick function meromorphic at x then g is analytic 
at x. This property propagates through a recursion to yield the n inequations 
fn+i{xj) 7^ Vj in Theorem [221 

Note 3.6 Reduction and augmentation at a point of analyticity are inverse oper- 
ations. 

More precisely, 

(1) if / G V is analytic at x and if g is the reduction of / at x then / is the 
augmentation of g at x by f(x), f'(x); 

(2) if g G P is analytic at x and if / is the augmentation of g at x by ao, ai then 
g is the reduction of / at x. 

However, when we admit augmentation at poles, it can happen that two different 
functions have the same augmentation at a point. The function z is the augmenta- 
tion of the zero function at by 0, 1, and is also the augmentation of the function 
-l/(2z) at by 0,2. 

Since reduction preserves V we may consider repeated reduction of a function 
/ G V. The process will terminate only if we reach a constant function. This can 
of course happen: if f(z) = l/(a — z) for some a G R and if x G M. \ {a} then 
/ G V and the reduction of / at x is the constant function a — x. 

4 An example 

A simple two-point interpolation problem will illustrate both the power of Nevan- 
linna's method and an important subtlety: 

Construct all functions f G V, analytic at 0, 1, such that 



If / is such a function then we may construct the reductions / 2 , f% of /, / 2 at 0, 1 
respectively. We find 



/(0) = 0, /(!) = !, f(0) = /'(!) = 2. 



(4.1) 



Hz) 



1 1 



m i 



f(z) 2z 



f(z) 2 2^2 
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so that / 2 (1) = -|, /^(l) = § and 



By Proposition 13.41 /■?. / 3 G "P. Furthermore, / 3 is analytic at 1 and, since /2(C)) 7^ 

00, 

/a(0) + -§, f 3 (l) ^ 00. (4.2) 

Conversely, suppose that / 3 G "P is meromorphic at 0, 1 and satisfies the in- 
equations (14. 2L We can augment fa, fa at 1,0 by (— 1,|),(0,2) to obtain f 2 , f 
respectively, and we obtain / G P satisfying the conditions (14. ip . We can write 
the general solution / of our problem as a continued fraction 

f (z) = -2 I 

2^ "r" 2 3(z-l) I 

parametrized by an arbitrary function f 3 G P, meromorphic at 0, 1 and satisfying 
the pair of inequations (14.21) . 

What if we choose f'3 G V that does not avoid the forbidden values - say 
/a(l) 7^ 00 but ^3(0) = — I? The augmentation process still produces an / G V, 
but it is not a solution of the interpolation problem, for 

/(0) = 0, /(1) = 1, /'(I) — 2 but/'(0)<2. 

The strict inequality arises because f 2 has a pole at 0, and so, according to Propo- 
sition [33(2), the augmentation / of f 2 at by 0, 2 satisfies /'(0) < 2. 

This example illustrates the need for the inequations f n +\(xj) 7^ yj in Theorem 
12.21 and so provides a counterexample to [TT], Satz I]. 

A second example illustrates another fact: if one performs Nevanlinna reduction 
on a solvable problem dNPV with distinct target values Wj one can obtain a 
problem in which some target values coincide, and consequently, after a second 
reduction, a problem in which 00 is a target value. 

Construct all functions f G V such that 

f(j) = fU)=j, J = 1,2,3. (4.3) 
After reduction at 1, we seek f 2 EV such that 

/ 2 (2) = / a (3) = 0, m) = i, m) = \. 

Since the target values at 2 and 3 are equal, the reduction / 3 of f 2 at 2 must satisfy: 
f 3 has a pole at 3 with residue —2. 
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5 Reduction and Schur complementation 



In this section we shall show that Nevanlinna's reduction process corresponds to 
Schur complementation of the corresponding Pick matrices. To this end we need 
to broaden the scope of the problem to encompass interpolation values at oo. 
Nevanlinna supposed in a footnote that the interpolation values Wj were distinct, 
but he ignored the possibility that equal target values might nevertheless arise in 
the course of the recursion (as they do in the example (14. 3ft above). 
Problem dNPVoo- Given distinct points X\, . . . ,x n G R and given w\, . . . , w n G 
RU {oo} and Vi, . . . , v n > such that Vj > whenever Wj = oo, determine whether 
there is a function f G V, meromorphic at x\, . . . , x n , such that, for j — 1, . . . , n, 

(1) f(xj) = Wj and f\xj) = Vj if Wj ^ oo, 

(2) / has a simple pole at Xj of residue —1/vj if Wj = oo. 

The Pick matrix of the data is the n x n matrix M = [my] where 

1 < i = j < n 
i 7^ j, Wi, Wj E R 

Wi = oo, Wj G R 

Wi G R, Wj = oo 
Wi = Wj = oo, i ^ j. 



rrti 



Vi if 

Wi - Wj 

if 
if 



X i 




Xj 




1 


X i 




Xj 




1 




Xj 




X i 










if 
if 



Note 5.1 Problems dNPV and dNPVoo are equivalent. 

That is, the solution sets of related problems are in bijective correspondence and 
the corresponding Pick matrices are diagonally congruent. 

For clearly every problem dNPV is also a dNPVoo, with the same Pick matrix. 
On the other hand, if / G V is a solution of Problem dNPVoo then h == — l/(f— a), 
for a G R \ {w\, . . . , w n }, is a solution of the problem dNPV with data 

x u . . . ,x n G R, wi(a), . . . , w n (a) G R, v x {a), . . . , v n (a) > 

where 



and 



«M)-{- w r a) *«y-i, (5 - i) 

= I Vj/(Wi ~ a? l WjER (5-2) 

{ Vj it Wj = oo. 

Conversely, if h G V satisfies the problem with these data then / = f a — 1/h solves 
the original problem dNPV'^. Furthermore, if the Pick matrices of the / data and 
the h data are M and M(a) respectively then 

M = DM(a)D, (5.3) 
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where D = diag{di, . . . , d n } and 

^ _ f w j ~ a if Wj GM. 

3 1 1 if U>j = oo. 

Thus the Pick matrices of the two problems are diagonally congruent. 
The following statement is simple to verify. 



(5.4) 
□ 



Proposition 5.2 Let g, h be solutions of problems dNPVoo with interpolation 
nodes x\, . . . ,x n , having corresponding Pick matrices M g ,Mh respectively. If g = 
c 2 h + d for some c, d G R with then 



M q = AM ft A 



where A = diag {c £l , . . . , c Sn } and 



1 ifg(xj) G 



' 3 I -i tfg(%j) = oo. 

We denote by schur A the Schur complement of the (1,1) entry of a matrix A. 

Theorem 5.3 Let Problem dNPVoo have Pick matrix M, and suppose that n > 1 
and Vi > 0. If f is a solution of Problem dNPVoo then the reduction g of f at Xi 
is a solution of the problem dNPVoo with data 

x 2 ,...,x n , w' 2 , ...,w' n G Ru{oo}, v' 2 ,...,v' n >0 



and Pick matrix M satisfying 



AM A = schur M, 



where A = diag {A2, . . . , A„} and 
The data w'j,v'j are given explicitly by 



W\ — Wj if Wi, Wj G R and w\ 7^ Wj 
1 otherwise. 



(5.5) 



(5.6) 



+ 



Wj — W\ V\(Xj — X\) 



w d = < 



vAxj — xi) 

1 



Vi(Xj — X\) 

00 



if Wj 7^ w\ in R 
if wi G R, Wj = 00 
if W\ = 00, Wj G R 



(5.7) 



V i 



1 



i/ Wj = wi, 
if Wj 7^ wi in R 



(wj — wi) 2 — £i) 2 

u,- ; ~ if Wi G R, = 00 or w\ = 00, iUj G 



if Wj = W1 . 



(5i 
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Proof. The appropriate values of wLv'; in equations (15.71) and (15. 8p are easily 
calculated from the relation 



1 

H 7 r II W\ E 



f(z)-w 1 Vi(z-Xi) 

f( z ) H 7 r if w i = °°- 

t>i(z - x x ) 



We first prove that AM A = schur M under the assumption that all Wj G R. Let 

M = K-]^ =1 , M = [m^ =2 , schur M = [m^. =2 . 
For j = 2, . . . , n, 



Uj 1 



m jj = V 'j = \ ( W j ~ ^l) 2 Vi(Xj — Xi) 2 



if lOj 7^ u>i 



if = Wx- 

If Wj = Wx then m^ = and \j = 1, while otherwise Aj = wx — Wj, and so 

, 2 ~ mxjirijx b 

A.m,-,- = m,-; — — = m- .. 

mn 

Thus AM A and schur M have the same diagonal entries. 

For 2<i^j<n, Wx, Wj ^ Wx, 

w'i- w 'i 1 / 1 1 1 1 

ra„ - 



Xi — Xj Xi — Xj \ Wi — Wx Vx(Xi — Xx) Wj—Wx Vx(Xj — Xx) 

b 



] / mi,-m,-i \ w,o- 







f my - 


m H mjx 






= 0,Ai 


= 1 and 



(wi - Wi)(wi - wj) V m n / h\ 

If ^ = Wx, Wj 7^ Wi then n 

1 m« 1 / m il m lj \ m 4 b . 

ray = = — = my = T-f-. 

Xi — xj Wi — Wj Wx — Wj \ mxx J \\j 

A similar calculation applies if Wi ^ wx,Wj = Wx, while if w.- t = Wj = wx then 

fhij = = mij = — 

AiAj 

Hence AM A = schur M when all Wj are finite. 

Now consider the general case of Problem dNPVoo, with some Wj possibly 
infinite. We reduce to the previous case by the transformation used in Note 15.11 
Choose a E R \ {w l} ■ ■ ■ ,w n } and consider the problem dNPV solved by h == 
— 1/(/ — a) at xx, ■ ■ ■ , x n \ this has data Xj, Wj(a) G R, Vj(a) given by equations 
(jBTTjl . ([E2D, and has Pick matrix M(ct) = D^MD" 1 with D as in (Q. Let /i 2 
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be the reduction of h at x±. If M(a) is the Pick matrix of h 2 at x 2 , • • • ,x n , then 
by the previous case, 

M( a ) = A'^schur M(a))A- 1 (5.9) 
where A a = diag {X 2 (a), . . . , A n (a)} and 



Xj(a) 



wi(a) — Wj(a) if wi(a) 7^ Wj(a) 
1 otherwise 

w i ~ w i t to / 

if lOi , Wj e K, Wi 7^ 



(wi — a)(tfj — a) 

if W! e R, wj = 00 

wi^— a (5.10) 

if W\ = 00, Wj e M 



Wj — a 



1 if Wi — Wj. 

Note that, if D 2 = diag {d 2 , ■ ■ ■ , <i n } then 

schur M(a) = schur (D -1 MD _1 ) = Dj^schur M)D 2 \ 
which, in conjunction with equation (15.91) . yields 

M(a) = A~ l D 2 l (schur M)D 2 1 A~ l . (5.11) 
A simple calculation gives the relation 

( h 2 1 

p = < [Wx - a) 2 Wx — a (5-12) 
I /i 2 + a if Wx — 00. 

Thus, in the case Wx — 00, 

M = M(a) = (iMa) -1 (schur M)(D 2 A a )- 1 . 

On comparing equations (15.4)) . (15.101) and (I5.6p we find that Z^A^ = A. Hence 
AMA = schur M when u>i = 00. 

In the case that wx £ K the functions g, /t 2 have poles at those Xj,j = 2, . . . , n, 
such that Wj — wx- By equation (15. 12)) and Proposition 15.2) 

M = A.V/(n)A 

where A = diag {(u>! — a) 62 , . . . , (u>i — a) £n } and 



-1 if Wj 7^ Wi 
1 if Wj = Wx- 



On combining this equation with (15. lip we obtain 

M = AA^D' 1 (schur M)D 2 ~ 1 A^ 1 A. 

By inspection, 

AA- 1 ^- 1 = A" 1 , 

and hence AMA = schur M, as required. □ 
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Note 5.4 For M^M as in Theorem \5.3[ 

~ _ ( t>r 1 det M ifwi = oo 

For 

det M = v\ det(schur M) = ^ det(AMA) = ^(det A) 2 det M = v x JJ A 2 (det M). 

3=2 

A result related to Theorem 15.31 was proved by I. V. Kovalishina [91 Section 
3]. She studied the analogous interpolation problem for several derivatives at a 
single point on the real axis. She gave a different formula for reduction, based on a 
"fundamental matrix inequality" , and showed that the Pick matrix of the reduced 
problem is congruent to the Schur complement of the (1,1) entry of the initial Pick 
matrix. Overall though there is not a great deal in common between our papers. 



6 Sarason's solvability criterion 

According to the Introduction of [12J, the most natural analogue of the stan- 
dard Nevanlinna-Pick problem for the version with boundary data is not Problem 
dNPV but rather the following "relaxed" version of it: 

Problem dNPV': As Problem dNPV, except that the condition f'(xj) = Vj is 
replaced by f'(xj) < Vj. 

We shall also need the relaxed version of the broadened problem, with infinities: 
Problem dNPV'^: As Problem dNPVoo except that the condition (1) f'(xj) = Vj 
is replaced by (l')f'(xj) < Vj if Wj G M., and condition (2) is replaced by 
(2') / has a simple pole at Xj of residue at most —1/vj if Wj = oo. 

Note 6.1 Problems dNPV and dNPV'^ are equivalent. 

For exactly as in Note IBTT| the transformation h = —l/(f — a), for suitable a G 
K, converts Problem dNPV'^ into dNPV , the two problems having diagonally 
congruent Pick matrices. 

Sarason proves his solvability result by first obtaining one for the relaxed prob- 
lem, a strategy that we shall follow. The following inductive proof is reminiscent 
of D. Marshall's elementary proof of Pick's theorem via Schur reduction |10j . 

Theorem 6.2 A problem dNPV or dNPV^ is solvable if and only if the corre- 
sponding Pick matrix is positive. 

The solution set of the problem, if non-empty, contains a real rational function 
of degree no greater than the rank of the Pick matrix. 

The problem is determinate if and only if the Pick matrix is positive and sin- 
gular. 
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Proof. In view of Note 16. 11 it suffices to consider Problem dNPV '. Let M be its 
Pick matrix. In the case that some v $ = the problem has a solution if and only 
if the constant function Wi is a solution, and it is easily seen that all statements of 
the theorem hold. We may therefore assume that all v% are non-zero. 

In the case n = 1, M = [v\] > and the problem has infinitely many solutions 

f(z) = w 1 + c(z - ari) 

for < c < vi, all real rational of degree at most 1. Thus the theorem holds when 
n = 1. 

Now consider n > 1, and suppose the theorem valid for n — 1. Again by Note 
16. 1^ the version for dNPV^ will also hold. 

Suppose that Problem dNPV has a solution /. The case in which the problem 
is solved by a constant function is easy, so we can assume that / is non-constant, 
whereupon all f'{xj) > 0. Consider the problem dNPV obtained when Vj is 
replaced by Uj = f'{xj) for j — 1, . . . , n. Then Uj < Vj and the Pick matrix M' of 
the modified problem satisfies 

M = M' + diag {v x - u x , . . . , v n - u n } > M' . (6.1) 

By Theorem 15. 3[ the reduction g of / at x\ is the solution of a problem dNPV^ 
onn-1 nodes having Pick matrix M satisfying 

M = A-^schur M')A _1 

for some nonsingular real diagonal matrix A. By the inductive hypothesis M > 0, 
and hence schur M' > 0. It follows from the identity ( 12. 3p that M' > 0, and hence, 
by the inequality (16. ip . that M > 0. Thus the existence of a solution of dNPV^ 
implies the positivity of the corresponding Pick matrix. 

Conversely, suppose that M > 0, and let r = rank M. Then schur M is positive 
and of rank r — 1. By Theorem 15.31 the problem dNPV^ with data 

x 2 ,...,x n , w' 2 , . . . ,w' n E RU {oo}, v f 2 ,...,v' n >0 

given by equations (15.7p . (15.8p has Pick matrix M congruent to schur M. Thus 
M > and rank M = r — 1. By the inductive hypothesis, the reduced problem 
dNPV'^ has a solution g G V that is real rational of degree at most r — 1. Let / 
be the augmentation of g at X\ by w%, V\ (it can be that x\ is a pole of g - see Note 
16.31 below). By Theorem 13.41 / G P and /(xi) = w%, f'(xi) < Vi, and by Note I3.3[ 
/ is real rational of degree at most r. We have 

J{Z)—W\ Vxiz — Xx) 

and hence, from condition (1), if Wj ^ W\ then 

1 _ 1 x , _ 1 

f(Xj) — Wi Vi(Xj — X\) 3 Wj — Wi ' 
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so that f(xj) = Wj, and 



P(xa) 1 ,. . 1 , Vj 



(Wj — Wi) 2 — Xi) 2 v l( x j — x i) 2 ( w j — w i) 2 

so that f'(xj) < Vj. If, however, j > 2, Wj = Wi and g has a simple pole of residue 
Rj < — 1 / -Uj at Xj then we have 

1 Rj ' Ofl) 



/(z) — Z — Xj 

for 2 in a punctured neighbourhood of X-,-, and hence f(xj) = Wj and 



f'( x j) = —5~ ^ 
Kj 

Thus / is a solution of the original problem dNPV'. Hence M > implies 
solvability of dNPV' by a real rational function of degree at most rank M. 

It remains to prove the statements about determinacy. Suppose M is singular 
and positive. If all the Wj are equal then the off-diagonal entries of M are all 
zero, and since M is singular some Vi is zero, and hence Problem dNPV has a 
unique, constant, solution. We may therefore suppose that wi / w 2 . Let fx, fi 
be solutions of Problem dNPV . Since M is singular, so is M, and hence by 
the inductive hypothesis the problem dNPV^ obtained by reduction at x\ has a 
unique solution g. Thus, for z£ll, 

1 1 . . 1 1 
+ 777- wT \ = <?(z) = -^73 — + 



h(z)-w x fx{xi){z - Xx) f 2 {z) -Wi ft(xi)(z - Xx)' 

Since W2 7^ w\ we may substitute z = X2 to obtain 

11 11 
+ ^7^7 r = + 



w 2 -wx fx(xx)(x 2 - xx) w 2 -wx f 2 (xx)(x 2 - Xx)' 

It follows in turn that f[(xx) = f 2 ( x i), that —l/(fx — Wx) = —l/(f 2 — wx) and that 
fi — /2- Thus if M is singular and positive then Problem dNPV' has a unique 
solution. 

Conversely, if M > then also M > and so, by the inductive hypothesis the 
reduced problem has two distinct solutions gx, g 2 6 V. The augmentations of gx, g 2 
at Xi by Wx,vx are distinct solutions of Problem dNPV. □ 

Note 6.3 Augmentation at poles and strict inequality in Problem dNPV . 

In the proof of sufficiency above we construct a function / e V by augmenting a 
rational function g by v 1 at the interpolation node Xx- In the event that x\ is a 
pole of g, according to Proposition 13. 1[ we have f'(xx) < Vx, with strict inequality. 
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We saw in Section H] that one can choose g to have a pole at x±. It can even be 
that g is unique and has a pole at X\. Consider the problem dNPV with data 



J 



1 2 3 

-1 1 

-1 1 

12 1. 



Here 



M 



1 1 1 
1 2 1 
1 1 1 



> 



rank M = 2. 



If we reduce at x 2 = we arrive at the problem of finding g 6 V such that 



The unique g G V satisfying these equations is g{z) = —l/(2z), which does have 
a pole at x 2 = 0. On augmenting g at by w 2 = 0, v 2 = 2 we obtain the unique 
solution f(z) = z of the problem dNPV , and as expected, 



Note 6.4 A statement of Nevanlinna. 

At what point of [H] did Nevanlinna err? Consider, in the light of the foregoing 
example, the statement [TU page 4]: "Notwendig und hinreichend dafiir, dafi eine 
Funktion der verlangten Art existiert, ist, dafi entweder die Gleichungen (7) beste- 
hen, in welchem Fall die funktion eine reelle Konstante ist, oder t\ = z[ > ist 
und eine Funktion f 2 (x) existiert, die der Bedingung (2) geniigt und in den (n — 1) 
Punkte (12) die Werte (13) mit den Ableitungswerten (14) annimmt." 

In the terminology of this paper, he is asserting that Problem dNPV is solvable 
if and only if either it is solved by a constant function or its reduction at an 
interpolation node is solvable. He evidently overlooked the possibility described in 
Note l6.3t in proving sufficiency one augments f 2 at xi, and x\ may be a pole of j%. 
In this case one obtains a strict inequality: f'{x\) < z[. Indeed, immediately before 
the above paragraph he wrote "Ferner bestatigt man unmittelbar, dafi die Werte 
der Funktion und ihrer Ableitung in den Punkten (3) mit den gegebenen Werten 
(4) bzw. (5) ubereinstimmen," a claim which fails in the case of augmentation at 
a pole. □ 

We now deduce Sarason's solvability theorem from Theorem 16.21 
Proof of Theorem 12.11 Necessity. Suppose that Problem dNPV has a solution 
/ £ V but that its Pick matrix M is neither positive definite nor minimally positive. 
A fortiori f solves Problem dNPV, and so, by Theorem 16.21 M > 0. Since M 
is not positive definite, M is singular, and so Problem dNPV' has the unique 
solution /. Since M is not minimally positive there is some index % and some 



0(±1)==F 



l 

2' 




1 

2' 



1 = f'( X2 ) <v 2 = 2. 
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positive Vi < Vi such that M' > 0, where M' is the matrix obtained when the (i, i) 
entry V{ of M is replaced by v/. Again by Theorem 16.21 there exists h G V such 
that h(xj) = Wj for each j, h'(xj) < Vj for j / i and h'(xi) < v{ < V{. In view 
of the last relation and the fact that f'(xi) = Vi we have h 7^ /, while clearly h 
is a solution of Problem dNPV' , as is /. This contradicts the uniqueness of the 
solution /. Hence if the problem is solvable then either M > or M is minimally 
positive. 

Sufficiency Suppose that M is minimally positive. By Theorem 16.21 there is an 
/ G V such that f(xj) = Wj and f'{xj) < Vj for j = 1, . . . , n. If in fact f'(xi) < Vi 
for some index % then consider the matrix M' obtained when the (i, %) entry V{ of 
M is replaced by f'{xj). Since M' is the Pick matrix of a problem dNPV' that is 
solvable (by /), we have M' > by Theorem 16.21 and so M majorises the non-zero 
positive diagonal matrix diag {0, . . . , v i — f'(xi), . . . , 0}, contrary to hypothesis. 
Thus f'{xj) = Vj for each j, that is, / is a solution of Problem dNPV. Moreover, 
since M is singular, / is the unique solution of Problem dNPV', hence is real 
rational of degree at most rank M. 

The above two arguments are essentially Sarason's. For the case M > we 
give a new, elementary, argument based on Nevanlinna reduction. We prove by 
induction a strengthening of the statement of the theorem. Let S(n) be the asser- 
tion: 

// the Pick matrix of Problem dNPV is positive definite, if T is a finite subset of 
R \ {xi, . . . ,x n } and y t G IR U {00} for every t G T then there are at least two 
functions f G V that solve Problem dNPV, are real rational of degree at most n 
and satisfy fit) 7^ y t for every t G T . 

The point of the condition f(t) 7^ y t is to ensure that all augmentations take 
place at points of analyticity. Once S(n) is established we can apply it with T the 
empty set to conclude the proof of the theorem. 

S(n) is true if n = 0. Here the Pick matrix condition is vacuously true, and one 
may choose any two distinct real-valued constant functions / with values outside 
the finite set {y t :tG T}. 

Let n > 1 and suppose that S(n — 1) is true. By the equivalence of Problems 
dNPV and dNPV oo the analogous statement will be true for dNPVoo. We apply 
reduction at x\. Consider the problem dNPV^ with data 

x 2 ,...,x n , w' 2 ,...,w' n , v' 2 ,...,v' n >0 

as in Theorem 15 .31 Let T 1 = T U {xi}, let y\ x = 00 and for t G T let 



y t -wx vxit-xx) 

with the natural interpretation y\ — 00 if y t — W\. By Theorem I5.3[ the Pick 
matrix M of this reduced problem is congruent to schur M, hence is positive 
definite. By the inductive hypothesis there are at least two real rational functions 
g G V of degree at most n — 1 that solve the reduced problem dNPVoo, are 
meromorphic at every point of T 1 and satisfy g(t) 7^ y\ for t G T 1 . In particular, 
g is analytic at x\. 
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For either g, let / be the augmentation of g at x\ by w\, v±: 

f(z) =w 1 + — . (6.2) 

Since g is analytic at X\, Theorem 13.41 and Note 13.31 tell us that / is analytic at x\, 
f(xi) = Wi, f'(xi) = V\ and / is real rational of degree at most n. For j — 2, . . . , n 
we have 

f 1 1 

, I + — if Wj + Wi 

g\Xj) = Wj = < Wj — w\ Vi(Xj — xi) 

^ oo if Wj = w\. 

If Wj 7^ W\ then the denominator in the right hand side of equation (16. 2p is non- 
zero at Xj and one sees that / is analytic at Xj and that f(xj) = Wj, f'(xj) = Vj. 
These relations also hold if Wj = w±. 

We claim that / is analytic on T and satisfies f(t)j^y t ,t&T. Indeed, since 
g is analytic at t and g(t) ^ yl, we have 

' -giz) 



Vi(z - x x ) 

is analytic and non-zero at t. Hence / is analytic at t. If y t ^ Wx, then since 
g{t) 7^ y} we have 

g(t) ^ + 77 7 

y t -wx vx(t-xx) 

and hence 

1 1 ,s , 1 



f(t)-Wx Vx(t-Xx) Vt-Wx 

and therefore fit) ^ y t . On the other hand, if it happens that w\ = y t , since g is 
analytic at t, it is clear that l/(/ — Wx) is analytic at t, hence fit) ^ w\. Thus 
f(t)^y t . 

In this way we construct at least two functions in V with the required properties. 
It follows by induction that S(n) holds for all non-negative integers n. □ 

There are other solvability criteria. Sarason gives a second one himself [T2| 
Theorem 2] based on an idea of Kotelyanskii, and D. R. Georgijevic Theorem 
1] gave another (for a more general problem) at about the same time. An earlier 
criterion may be found in [6l Theorem XIII. 1]. However, Sarason's first criterion, 
as in Theorem 12.11 is the simplest and most elegant. 

Note 6.5 Solvability of Problem ON PVoq. 

The above proof, with minimal changes, shows that Theorem 12.11 remains true if 
dNPV is replaced by dNPV^. 

Note 6.6 Functions not analytic at the interpolation nodes. 

Sarason and other authors (e.g. [3]) study problems like dNPV with less stringent 
regularity assumptions than analyticity at the interpolation nodes. Let us say that 
a function / e V is a weak solution of Problem dNPV' if, for 1 < j < n, f has 
non-tangential limit f(xj) = Wj and angular derivative f'(xj) at most Vj at Xj. 
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Proposition 6.7 The following five statements are equivalent: 

(1) Problem dNPV' has a weak solution; 

(2) Problem dNPV has a solution; 

(3) Problem dNPV' has a rational solution; 

(4) Problem dNPV' has a real rational solution; 

(5) The Pick matrix of the problem is positive. 



Proof. By TheoremES (2)^ (3) <£> (4) (5), and obviously (2)^(1). It is also 
true that (1)=^(5). For suppose that / is a weak solution of Problem dNPV . It 
is simple to show that /' has non-tangential limit equal to the angular derivative 
f'(xj) at Xj (for example, [121 VI-1]). For any e > let f £ (z) = f(z + ie). Clearly 
f £ EV and f £ is analytic on JR.. By Theorem 16.21 the Pick matrix M(e) = [m,ij(e)] 
is positive, where 

mijie) = { f e (xj) - f e (xj) . f . ^ . 

Ob j Ob q 



Since 



3 



lim f' e (xj) = lim f'(xj + ie) = f(xj) < Vj, 



and f e (xj) — * vjj as e — > 0, it follows that M(e) tends to a matrix M' > 0, where 
the Pick matrix M of Problem dNPV' satisfies 

M = M' + diag { Vl - f'(x 1 ), ...,v n - f'(x n )}. 

Hence M > 0. Thus (1)^(4). □ 
In view of these equivalences there is no loss (as far as existence goes) in restrict- 
ing ourselves to interpolation by functions that are analytic at the interpolation 
nodes. 



7 Parametrization of all solutions 

We turn to the proof of Theorem 12.21 The method of proof is essentially that of 
Nevanlinna, but in adding the necessary considerations of infinities and augmen- 
tation at poles one cannot avoid a modicum of algebraic detail. 

Lemma 7.1 Let the problem dNPVoo with data 

x 1 ,...,x n eR, wx, . . . ,w n e KU {oo}, vi,. . . ,v n > 
have positive definite Pick matrix. Define 

<...,<6MU{oo}, >0, fork = l,...,n, 
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recursively by the relations 



Wj = Wj, 



v 3 = v h 



forj = l,...,n 



and, for 2<k + l<j<n, 



(7.1) 



w 



fe+i 



+ 



w k_ w k v%(Xj—Xk) 



,fc+l 



3 V k( X 3 - X k) 

oo 

1 



if wjel and Wj = oo 

if w*, w'eR and Wj ^ w\ 

if w% = oo and w'gI 

if w \ = w ^ 



(7.2) 



if G M, = oo or u>|? = oo, Wj G 



if w'el and 7^ w\ 
if wl = w^. 



Suppose further that y\ , . . . , y%-i £ SU are defined for k = 2, 
Dk-i = 00 an d, for 1 < j < k < n, 



(7.3) 
n + 1 by 



+ 



1 



Vj-K v^Xj-Xk) 



if wl G 



1 



(7.4) 



For k — 1, . . . ,n, if a function f G P satisfies 
(1) f is a solution of problem dNPVoo with data 
x k ,...,x n eR, ^,...,^elU{oo}, 



«* ...,^>0, (7.5) 



(2) / is meromorphic at x lt . . . , x^-i and /(x,) 7^ n^ /or j — 1, . . . , k — 1 
i/ien £ne reduction g of f at Xk satisfies 
(!') g is a solution of problem dNPVoo with data 



•Ek+li • • • 1 %n G 



fc+1 



fc+1 , . . . , w n 



wt +1 G RU {00}, 



,fc+i 



V k+1 1 ■ ■ ■ i V n 



fc+1 



>0, 



(2') g is meromorphic at Xi, . . . , x^ and g(xj) 7^ y- /or j = 1, . . . , k. 

Conversely, if g G V satisfies conditions (V) and (2') then the augmentation f of 
g at Xk by w^.,v^. satisfies conditions (1) and (2). 
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The expressions on the right hand side of equation (17.41) have their natural inter- 
pretations in the event that y k = oo or y k = w*. 

Proof. The proof is a re-run of the proof by induction of S(n) in the final part 
of the proof of Theorem 12.11 but with extra detail. If is the Pick matrix of the 
interpolation problem with data (17.51) . k = 1, . . . , n, then for k < n — 1, Mk+\ is 
congruent to schur M^. Since M\ > it follows that Mk > for k — 1, . . . , n. In 
particular v k > for 1 < k < j < 



n. 



Fix fc, 1 < k < n. Suppose that / G V satisfies conditions (1) and (2) and that 
g is the reduction of / at Xk- By Theorem 13.41 g G V and g is analytic at Xk- We 
shall show that g satisfies (1') and (2'). 

Consider first the case that w k = oo. Since / has a pole at Xk with residue 
— l/v k we have 

v k {z - x k ) 

Clearly g is meromorphic at x\, . . . ,x n . Consider g(xj) for j = k + 1, . . . , n. If 
Wj G M then f(xj) = w k and so 

while if w k = oo then / has a simple pole at Xj with residue —l/v k , and hence the 
same is true for g. Since v k+l = v k when w\ = w k = oo it follows that g has a 
simple pole with residue — l/v k+1 at Xj. Hence g satisfies (1') when w k = oo. 

Since g is analytic at Xk we have g(xk) ^ oo = Vk +1 - For j — 1, . . . , k — 1 we 
have /(x-,) 7^ and so 



1 nixj-xk) 1 

(note that this assertion remains valid when y k = 00). Thus g satisfies (1') and 
(2') when w^ = 00. 

Now consider the case that w^ G M. We have 

0(*) = -7TT — + ^ V- ( 7 - 6 ) 

f{z) - W% Vl{z - Xk) 

For j — k + 1, ... ,n we consider the subcases (i) w k G K, 7^ u^, (ii) wj 1 = G 
K, (iii) = 00. In subcase (i) we have 



g(Xj) = r — r + -r- r = W k+l , 

Wj - W% V%{Xj - x k ) 



g\xj) 



v ) 1 _ M 



V : 



(w k - w k ) 2 vl(xj-Xk) 2 3 



and so g satisfies (1') at Xj. In subcase (ii) one sees from equation (17.61) that g has 
a simple pole at Xj with residue —l/v k . Since in this subcase w k+1 = 00, v k+1 = v k 
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it is again true that (1') holds at Xj. In subcase (iii) / has a pole of residue — l/uj 



at Xj, and we find that 



— u k+1 

V%{Xj-X k ) J ■ ~ ' " J V%[Xj-X k f 



Hence g satisfies (1') when w\ G R. 

We must show that g satisfies (2') when G R. Again g is meromorphic at 
x\,...,x k and g{x k ) 7^ oo = y k +1 - For j = 1, . . . , k — 1 we have by hypothesis 
/(xj-) 7^ j/J. If y) ^ w£ then 

f ( x i) ^ — k + -TT~ V = 

while if = w£ then <?(xj) 7^ 00 = y k j +1 ■ Thus g satisfies (2') as required. We 
have shown that in all cases g satisfies both (1') and (2'). 

Conversely, suppose that g G V satisfies conditions (1') and (2'), and let / be 
the augmentation of g at x k by wj^,v k . We shall prove that / satisfies conditions 
(1) and (2). 

Consider the case that w k = 00. We have 

/(*) = g{z) 



v%{z-x k y 



Clearly / is meromorphic at X\, . . . ,x n . Since g(x k ) 7^ y^ +1 = 00, g is analytic at 
x k . Hence / has a pole at x k with residue —l/v k , so that / satisfies (1) at x k . For 
j — k + 1, . . . , n we have 



g(x j )=wf+ 1 . (7.7) 

Either (i) w k+l = 00 or (ii) w k+1 G R. Case (i) occurs when w k = Wj. Here /, g 
both have poles at Xj with residue — l/vj +1 . Since v k+1 = v k when wfi = w k we 
deduce that / has a pole at Xj with residue —1/vf-. In Case (ii) 

f( Xj ) = r = w ;el, 

fc+i j _ „, fc 



- ^fc) 2 



7 ' 



Thus / satisfies condition (1) when w k = 00. 

For j = 1, . . . , k-1 we have g-(^) 7^ ^ fc+1 . Either (i) = 00 or (ii) G R. 
In Case (i) g is analytic at Xj, and hence / is also, whence f(xj) 7^ 00 = y k -. In 
Case (ii) 

f(xj) * y- +1 - ■ — 7 = yf • 

Thus / satisfies (1) and (2) when u>£ = 00. 
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The remaining possibility is that w k G M. Here 

1 1 



f(z)-w k v k (z-x k ) 



g(z). (7.1 



Since g is analytic at x k we have f(x k ) = w k , f'{x k ) = v k by Theorem 13.41 For 
j — k + 1, . . . , n either (i) w k+l = oo or (ii) w k+1 G K. In Case (i) g has a pole at 
Xj with residue — l/v k+1 , and equation (l7.8p shows that f(xj) = w k , f'(xj) = v k+1 . 
It is clear from equations (17.21) and (17. 3p that w k = w k and v k = Vj, and so 

ffa) = w k , fix,) = v). 
In Case (ii) one of the first two alternatives in equation (17.21) holds, and so either 

w k = oo w k+1 = v k+1 = v k 

3 3 vl(Xj-X k y 3 V^Xj-Xk) 2 



or 



w k el,^ / w k , w) 



fc+i _ 1 1 ^fc+i 



w k - W k V%(Xj - x k ) ' 3 {w k - W^) 2 V^(Xj - x k ) 2 



In the former case it is clear from equation (17. 8p that l/(/ — w k ) vanishes and has 
derivative 



' v k+1 = -v k 



V k (Xj - x k ) 2 3 3 

at Xj, which implies that / has a pole at Xj of residue — l/v k . 
In the latter case (w k 6 t, ^ w k ), 



w k+1 



f(Xj) - W k Vk( X j ~ X k) 3 Wj - w k 

and so f(xj) = w k . On differentiating equation (17. 8p and setting z = Xj we have 



ffa) 1 V 3 



k 



(w k - u>£ ) 2 v\(xj-x k ) 2 3 (Wj - w k ) 2 

and so f'(xj) = v k . Thus in all cases / satisfies (1). 

We must show that / satisfies (2) when w k G R. For j = 1, . . . , k — 1 we have 

9 fa) Vj +1 and 

ii i 

fc+i 



77"^ k = — f 7 - 9 fa) ? ^ - W - ( 7 - 9 ) 

ffa)-w k V k ( Xj -X k ) V k (Xj-X k ) 3 

If y k+1 = oo then y k = w k . Since l/(/ — w k ) ^ oo at Xj we have ffa) ^ = y k . 
On the other hand, if y k+l G M then equation (17. 9p becomes 

1 + 1 



ffa)-™* v*j- wk k 
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from which it follows that f(xj) 7^ y*. This concludes the proof that / satisfies 
conditions (1) and (2). □ 
Note that if one starts the recursion of Theorem 12.21 with a function f n+ \ G V 
that is meromorphic at each Xj but does not satisfy the inequations f n+ x(xj) 7^ yj 
then one obtains a function / G V that is a solution of Problem dNPV' , not 
dNPV. 

We can give an alternative expression for the parametrization of solutions in 
Theorem 12.21 in terms of a linear fractional transformation (as in [TH §4] , where 
Nevanlinna considers the generic case, in which all s k are finite). For any 2x2 
matrix A = [a^] let us denote the corresponding linear fractional transformation 
by L[A]: 

L\A \{w) — , wGLU 00 . 

1 n ' a 21 w + a 22 1 J 

In this notation the relationship between f k+ i and its augmentation f k at x k by 
s k ,t k can be written 

f k (z) = L[A k (z)](f k+1 (z)) 



where, for k 



Adz) 



t k (z-x k ) 1 

t k (z-x k ) 



S k t k (z — x k) —tk( z ~Xk) — S k 

t k (z-x k ) -1 



if s k 



00 



(7.10) 



if s k G 



Note that, in either case, det A k (z) = t k (z — x k ) 2 . The recursion for / in Theorem 
12.21 becomes: 

f(z) = f l {z) = L[A 1 {z)}{f 2 {z)) = L[A l {z)A 2 {z)}{h{z)) = ... 
= L[A 1 (z) . . . A n (z)](f n+ i(z)). 

We therefore arrive at the following linear fractional parametrization. 

Corollary 7.2 If Problem dNPV has a positive definite Pick matrix then its gen- 
eral solution is 

= a(z)h(z) + b(z) 
J{ ' c{z)h{z) + d(z) 
where a, b, c, d are real polynomials of degree at most n satisfying, for some m > 0, 



(ad — bc)(z) = m Y\( z ~~ x kY 



k=l 



A 1 (z)A 2 (z) ...A n (z), 



and given by 

a(z) b(z) 
c(z) d(z) 

where A k (z) is given by equations (17.101) . the quantities s k ,t k and y k are as in 
Theorem 12.21 and h is an arbitrary function in the Pick class that is meromorphic 
at Xi, . . . , x n and satisfies h(xj) 7^ yj, j = 1, . . . , n. 
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Yet another way to write down the parametrization of solutions in Theorem 12.21 
IS cLS cL continued fraction. In the event that si, . . . , s n are all finite we have, for 
1 < k < n, 

fk{z) = s k + 



z^r - fk+l( Z ) 



and so 



t k (z-x k ) 



f(z) = f 1 (z) = s 1 + ' 



Si + 



17(^7) _ ^ 
1 1 



Continuing in this fashion we obtain the following statement. 

Corollary 7.3 Let Problem dNPV have a positive definite Pick matrix and let 
Sk, tk andi/k be defined as in Theorem \2.2\ for k = 1, . . . , n. Suppose that s\, . . . , s n £ 
R. The general solution of the problem is 

f( z ) = Sl + 1 ~ i • • • i h(z) 

where h is an arbitrary function in the Pick class that is meromorphic at x±, . . . , x n 
and satisfies h{xj) 7^ yj, j = 1, . . . , n. 

It is of course simple in principle to modify this expression for the case that some 
Sk are infinite. 

We note a relation between the quantities Sj,tj etcetera that appear in the 
parametrization theorems. If one iterates the determinantal formula in Note 15.41 
one obtains 

ht 2 ...t n Y[ (wj - s k ) 2 = det M 

l<k<j<n 

where the prime on the product sign indicates that any factor that is infinite or 
zero is to be omitted. 



8 Boundary interpolation in the Schur class 

A version of Sarason's Theorem also holds for boundary interpolation by functions 
in the Schur class. Its derivation from Theorem 12.11 with the aid of Cayley trans- 
forms, is straightforward and is outlined by Sarason [12], so we shall merely state 
the result and go on to describe the parametrization of all solutions, which is easily 
derived from Theorem 12.21 We denote by T the unit circle in the complex plane. 
Problem dNPS: Given distinct points £1, . . . , £ n G T and given rji . . . , r) n £ T and 
Pi-, ■ ■ ■ 1 Pn > 0, determine whether there is a function if £ S, analytic at £1, ... , £ n; 
satisfying 

<P(t,j) = Vj, OVXf j) = Pj, j = 1, . . . , n, (8.1) 
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and if there is, describe the class of all such functions <p>. 

The Pick matrix of the problem is the matrix M = [mjj]Jj_i where 

f 1 - ViVj r ■ , ■ 
=- jL % J % T J> 
1 - iiii 
Pj if i = 3- 

Here A(p(£) denotes (for a function (p different iable and non-zero at £ = e 1 ' G T) 
the rate of change of the argument of (^(e 1 *) with respect to t, so that for (eT 

In particular, if <p G S and |<£>(0I = 1 at some £ G T then ^(e 1 *)) 2 has a local 
maximum, hence a critical point at £, from which it follows that £</?'(£)v(0 ^ ^- 
Hence, for such <£> and £ we have 

{Am) = W' (8,2) 



Theorem 8.1 Problem dNPS is solvable if and only if the corresponding Pick 
matrix is positive definite or minimally positive. 

The solution set of the problem, if non-empty, contains a finite Blaschke product 
of degree no greater than the rank of the Pick matrix. 

The problem is determinate if and only if the Pick matrix is minimally positive. 

For r G T we denote by C T the Cayley transform 

aw = ii±± 

T — A 

C T maps B to II and T to R U {oo}. The following is a simple calculation. 

Lemma 8.2 Let (p G <S, let r, a, £, 77 G T and suppose that 

T£ = e 2 ^ ^ 1, ar/ = e 2ie ^ 1. 

Suppose that (p is analytic at £ and satisfies (p(£) — r\. If f = C a o ip o C7 1 i/ien 
f <EV, f is analytic at C T (£) and 

^(o = ^r°a(o- 

sm -0 

Theorem 8.3 Suppose that the Pick matrix of Problem dNPS is positive definite. 
Let r,a ET be distinct from £ 1; . . . , £ n and f]x, ■ ■ ■ ,f] n respectively. There exist 

Si,...,s n 6lU{oo}, >0, 7/1, . . . ,y n G MU {00} 

swc/i that the general solution of Problem dNPS is <p = C~ l o f\ o C T where the 
functions f n +i, ■ ■ ■ , fi G V are given recursively by: 
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(1) f n+ \ is any function in V that is meromorphic at C T (£i), . . . , C T (£ n ) and 
satisfies f n+1 o C T (Q ^ y jt j = 1, . . . , n; 

(2) fk is the augmentation of fk+i at C r (£fc) by s&, /or fc = n, . . . , 1. 

Proof. With the aid of Lemma 18.21 one sees that a function ip G 5 is a solution 
of Problem dNPS if and only if C CT o o C~ x is a solution of Problem dNPV with 
data 

where, if f Q = e 2l ^ j and cx?7j = e 2l6j , 

sin 2 ■0,- 

v 3 = ■ 2 n Pi> j = 1, . . . ,n. (8.3) 
sin &j 

By Theorem 12 . 21 there exist Sj G MU{oo}, ^ > 0, yj G MU{oo} for j = 1, . . . , n such 
that the recursion (l)-(2) parametrizes all solutions of this problem dNPV. We 
then obtain the general solution ip of Problem dNPS by taking C a oipo = f\. 
□ 

Note that we can calculate the Sj,tj and yj explicitly by substituting Xj = 
C T (£j), Wj = C a (r)j) and Vj as given by equation ( 18. 31) in the formulae in Lemma 
17.11 We can also interpret the result as a linear fractional or continued fraction 
parametrization of all solutions, as in Corollaries 17.21 and 17.31 

Note 8.4 Reduction and augmentation in the Schur class 

It might seem more natural to parametrize the solution set of Problem dNPS by 
elements of S satisfying suitable inequations, rather than elements of V. This is 
indeed possible, but reduction and augmentation, when transferred to the Schur 
class by Cayley transformation, are very cumbersome. We therefore focus in this 
paper on the Pick class, even though our original goal was to understand problem 
dNPS and despite the modest complications required to deal with possible infinite 
target values (which do not arise in the case of S). Similar considerations apply 
to boundary interpolation problems in which values of higher derivatives of the 
function are prescribed. We intend to show in a future paper how Nevanlinna 
reduction affords an elementary solution of such problems as well. 
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